First-principles study of pressure-induced phase transition and electronic property of 
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We have performed a systematic first-principles investigation to calculate the structural, elec- 
tronic, and magnetic properties of PbCr03, CrPbC>3 as well as their equiproportional combination. 
The local density approximation (LDA)+f/ and the generalized gradient approximation-!-?/ theo- 
retical formalisms have been used to account for the strong on-site Coulomb repulsion among the 
localized Cr 3d electrons. By choosing the Hubbard U parameter around 4 eV within LDA+?7 ap- 
proach, ferromagnetic, and/or antiferromagnetic ground states can be achieved and our calculated 
volumes, bulk moduli, and equation of states for PCO-CPO in R3 phase and PCO in Pm3m or R3c 
phases are in good agreement with recent experimental Phase I and Phase II [W. Xiao et at, PNAS 
107, 14026 (2010)], respectively. Under pressure, phase transitions of R3 PCO-CPO to PmZm PCO 
at 1.5 GPa and i?3 PCO-CPO to i?3c PCO at -6.7 GPa have been predicted. The abnormally large 
volume and compressibility of Phase I is due to the content of CrPbOs in the experimental sample 
and the transition of PbOe/2 octahedron to CrOg/2 upon compression. Our electronic structure 
study showed that there will occur an insulator-metal transition upon the phase transitions. Clear 
hybridization of Cr 3d and O 2p orbitals in wide energy range has been observed. 

PACS numbers: 61.50.Ah, 71.15.Mb, 75.25. +z 



I. INTRODUCTION 



Strongly correlated electron systems of transition- 
metal oxides with ABO3 cubic perovskite or pseudo cu- 
bic perovskite structures exhibit particular interesting 
physical properties. [3-0] Their ferroelectric, ferromag- 
netic, ferroelastic, multiferroic, and/or magnetoresistivc 
features originate from the mutual interplay of various 
degrees of freedom, including lattice, spin, charge and 
orbital, in their partially filled B site 3d electrons. The 
multiple chemical characters of the A ion with lone pair 
electrons, especially for Bi 3+ and Pb 2+ , also play an im- 
portant role. Hill Correctly describing of their electronic 
and magnetic structures are critical. 



In late 1960s, a few groups [7H12| successfully syn- 
thesized some perovskites containing Cr 4+ (CaCrOs, 
SrCrOa, and PbCrOa) under high temperature and high 
pressure. For PbCrOa (PCO), a lattice constant of about 
4.00 A for the cubic structure was determined by X-ray 
diffraction on single crystal and powder samples and pow- 
der neutron diffraction @-H[. It was reported that the 
PCO is an antiferromagnetic (AFM) G-type semiconduc- 
tor with a magnetic moment of ~1.9 fis on each Cr ion 
@, |H and with 0.27 eV activation energy. @ The Neel 
temperature (Tn) of about 240 K was obtained by exam- 
ining the magnetic and transport properties. 0, Q For 
CaCr03 and SrCr03, preliminary structural, magnetic, 
and conductive properties were also investigated. 
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12] From then on, in a long period of more than thirty 
years little works had focused on these systems due to 
difficulty of synthesis. However, a renewed interest on 
them has been inspired on their transport and magnetic 
properties, insulator-metal transition, and pressure be- 
haviors. [13-QjI Anomalous properties of Seebeck coeffi- 
cient, thermal conductivity, magnetic susceptibility, and 
room-temperature compressibility have been observed for 
SrCrOs by Zhou et al. They concluded that SrCrOs 
is nonmagnetic (NM) insulator and CaCrOs is also an 
insulator at low temperature. But more recent stud- 
ies claimed that both these perovskites are AFM met- 
als. [l8l [n| An orbital ordering transition from t\ g to 
dxyidxzdyz) 1 and electronic phase coexistence of C-AFM 
tetragonal and NM cubic phases have been discovered in 
SrCrOs. [lH Komarek et al. reported that CaCrOs is 
an intermediately correlated metal with similar C-type 
AFM ground state. Thus, whether these systems are 
metallic, stron gly correlated, and spin ordered is still con- 
troversial. [2a|2i! 

Recent experimental works of PCO concentrated on 
its structure, electron energy loss spectroscopy, mag- 
netic structure, and high pressure phase transition. 

Electron diffraction and high-resolution electron 
microscopy study revealed that the microstructure of 
"PbCrOs" is a rather complex perovskite with a compo- 
sitionally modulated a p x 3a p x (14~18)a p superlattice 
structure, where a p = 4.002 A is the lattice constant of 
the cubic PCO perovskite. [HI The magnetic structure 
of PCO is also complex. Alario-Franco et al. reported 
AFM ordering of the chromium moments at Tn~245 K 
with a spin-reorientation at temperature range of 185 K 
to 62 K and their magnetic hysteresis loops for "PbCrOa" 
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suggested weak ferromagnetism at low temperature. [151 ] 
Their resistivity measurements indicated two activation 
energies ranges with 0.11 and 0.26 eV in different tem- 
peratures. As for pressure study, recent work performed 
by Xiao et al. [16| observed a large volume collapse in 
the isostructural transition of cubic PCO perovskite at 
~1.6 GPa from Phase I to Phase II. They concluded that 
the transition seems not related with any change of elec- 
tronic state, but probably has tight relation with the ab- 
normally large volume and compressibility of the Phase I. 
The real Phase I might be a kind of mixture of PbCrC>3- 
CrPb0 3 (PCO-CPO) combination due to the fact that 
the cubic lattice constant is enlarged if the Cr0 6 / 2 octa- 
hedron could be replaced by Pb0 6 / 2 16}. 

In present study, we focus our sight on PCO, PCO- 
CPO, and CrPb03 (CPO) in the cubic perovskite struc- 
ture (space group PmZra) and some possible distorted 
perovskite structures, such as i?3c, i?3, and Pi/mmm 
phases. Electronic and magnetic properties as well as 
pressure behaviors have been systematically investigated 
by the first-principles electronic structure calculations 
based on density functional theory (DFT) and DFT+t/ 
schemes due to Dudarev et al. [22j The validity of the 
ground-state calculation is carefully tested. Our cal- 
culated lattice parameter and bulk modulus B for cu- 
bic PCO are well consistent with previous local den- 
sity approximation (LDA) and generalized gradient ap- 
proximation (CCA) results [ 1 ( { . The total energy, lat- 
tice constant, bulk modulus B, and spin moment of 
Cr ion for NM, ferromagnetic (FM), and AFM phases 
calculated in wide range of effective Hubbard U pa- 
rameter are presented and our calculated results within 
LDA+ U with [7=3-4 eV for PCO in Pm3m or i?3c 
phases and PCO-CPO in R3 phase accord well with ex- 
perimental [16J Phase II and Phase I, respectively. Our 
calculated spin moment by LDA+J7 is in good agree- 
ment with recent experimental value of saturation mo- 
ment M sa t=1.70 /J.B, which is deduced from the effec- 
tive moment M e ff=2.51 [l5[ according to the relation 
A/ off =2[5(5 + l)] 1 /2 = [M sat (M sat +2)] 1 /2. The P-V rela- 
tions of PCO and PCO-CPO are calculated to compare 
with experiment. The insulting property of "PbCrOa" 
at ambient condition is successfully predicted. 



II. COMPUTATIONAL METHODS 

First-principles DFT calculations on the basis of the 
frozen-core projected augmented wave (PAW) method 
of Bfochl [231 ] are performed within the Vienna ab ini- 
tio simulation package (VASP) [24|, where the exchange 
and correlation effects are described by the DFT within 
LDA and GGA [25|,|26j. For the plane- wave set, a cutoff 
energy of 500 eV is used. The fc-point meshes in the 
full wedge of the Brillouin zone (BZ) are sampled by 
12 x 12 x 12 and 6x6x6 grids according to the Monkhorst- 
Pack [H] scheme for PCO, CPO, and PCO-CPO in their 
cubic and rhombohedral unit cells, respectively. The cu- 



bic perovskite structure are built for nonmagnetic (NM) 
and ferromagnetic (FM) calculations and the rhombo- 
hedral unit cells (Fig. 1) for G-type antiferromagnetic 
(AFM) configuration, which is (4, i, |) order in terms 
of the original perovskite cell. The rhombohedral unit 
cell for PCO in its cubic structure is constructed using 
the i?3c space group (a=60°) with Pb atoms in 2a(0, 
0, 0) site, Cr in 2a(±, ±, ±) site, and O in 66(±, 0, §) 
site. Moving away slightly atoms from these high symme- 
try positions will lower the symmetry and consequently 
build the cell in real R3c or R3 phase. In this study, the 
Pb 5d 10 6s 2 6p 2 , Cr 3d 5 4s\ and O 2s 2 2p 4 orbitals are ex- 
plicitly included as valence electrons. The strong on-site 
Coulomb repulsion among the localized Cr 3d electrons 
is described by using the formalism formulated by Du- 
darev et al. [22|. In this scheme, the total LDA (GGA) 
energy functional is of the form 

-ElDA(GGA)+U = -E-LDA(GGA) 

a 

where p a is the density matrix of d states with spin 
cr, while U and J are the spherically averaged screened 
Coulomb energy and the exchange energy, respectively. 
In this work, the Coulomb U is treated as one variable, 
while the parameter J is set to 0.5 eV. Since only the 
difference between U and J is meaningful in Dudarev's 
approach, therefore, we label them as one single param- 
eter U for simplicity. To obtain the energy data under 
different pressures, we perform the structure-relaxation 
calculations at a series of fixed volumes. The correspond- 
ing pressure values are deduced from the energy-volume 
data by P=-dE/dV. 

III. RESULTS 

A. Structural and magnetic properties of cubic 
PCO 

Both spin-unpolarized and spin-polarized calculations 
are performed for cubic PCO. For NM, FM, and AFM 
configurations, the total energies (-39.685 eV, -39.706 
eV, and -39.818 eV, respectively) calculated within the 
DFT (t/=0) have no visible differences. After turning 
on the Hubbard U, the NM phase is not energetically fa- 
vorable both in the LDA+ U and GGA+ U formalisms 
compared with FM and AFM phases. In the follow- 
ing, we only present results of FM and AFM configura- 
tions. The dependences of total-energy differences (per 
formula unit at respective optimum geometries) between 
FM and AFM (£"fm — -Eafm), lattice parameter, bulk 
modulus, and spin moments of Cr ions on U for cubic 
PCO are shown in Fig. 2. The theoretical equilibrium 
volume, bulk modulus B, and pressure derivative of the 
bulk modulus B' are obtained by fitting the third-order 
Birch-Murnaghan equation of state (EOS) [28j. For com- 
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parison, recent experimental values [16| of ao and B for 
Phase I and Phase II as well as the experimental results 
of spin magnetic moment of Cr ions are also shown. In 
LDA or GGA (U=0 eV), the total energy of AFM phase 
is lower than that of FM phase. However, as shown in 
Fig. 2(a), the total energy of the FM phase decreases to 
become lower than that of the AFM phase when increas- 
ing U. At a typical value of U =4 eV, the total-energy 
differences between FM and AFM (Efm—Eafm) reach 
their minimums of —74 meV and —91 meV within the 
LDA+ U and GGA+ U formalisms, respectively. Note 
that we have also considered other types of AFM config- 
urations and our results show that the G-AFM in (i, |, 
^) order is the most stable state. These results are con- 
sistent with recent experimental observations [l5| , where 
they reported the AFM ordering of the chromium mo- 
ments at Tn^245 K with a spin-reorientation at temper- 
ature range of 185 K to 62 K and their magnetic hystere- 
sis loops for "PbCrCV suggested weak ferromagnetism 
at low temperature. As clearly shown in Fig. 2(b), our 
calculated lattice parameters of cubic PCO within two 
DFT+ U approaches are all by lar ge s maller than the 
experimental value ao=4.013 A pzl— Tol Il6| of Phase I. Cor- 
responding results of bulk modulus are dramatically big- 
ger than that obtained by experiment (B=59 GPa) for 
Phase I [Fig. 2(c)]. But the experimental lattice param- 
eter ao=3.862 A and bulk modulus 5=187 GPa of Phase 
II well lies in the range of our calculated results of cu- 
bic PCO, which supports their conclusion of cubic PCO 
perovskite [l6| that the real Phase I might be a kind of 
mixture of random PCO-CPO combination. In next sub- 
section, we will discuss carefully the structural properties 
of PCO-CPO. 

As shown in Fig. 2, the tendencies of ao, B, and spin 
moments of Cr ions for FM phase with U are similar 
to that of the AFM phase. Results of NM phase (not 
presented) further indicates its unstable nature at low 
temperature. As shown in Fig. 2(b), LDA underesti- 
mates the lattice parameter with respect to the experi- 
mental value while GGA coincides well with experiment. 
Our LDA and GGA results for NM phase accords well 
with previous first-principle calculations [161 ]. Result of 
LDA is due to its typical overbinding character. The 
LDA+ U method will lead to relative larger equilibrium 
volume compared to the LDA and therefore improves the 
agreement with experiment, especially for FM and AFM 
phases. At a typical value U =4 eV, the LDA+ U gives 
a =3.822 A for AFM cubic PCO which is very close to 
the experimental value. On the other hand, the GGA+ U 
enlarges the underbinding effect with increasing Hubbard 
U. As a comparison, at U =4 eV, the GGA+ U gives 
ao=3.921 A. As for the dependence of bulk modulus B 
on U shown in Fig. 2(c), it is clear that the LDA re- 
sults (230-153 GPa) are always higher than the GGA 
results (187-113 GPa), which is due to the above men- 
tioned underbinding effect of the GGA approach. At a 
typical value t/=4 eV, the LDA+ U and GGA+ U give 
5=198 (#=4.6) and 148 GPa (#=5.1) for AFM cubic 



PCO, respectively. Obviously, result calculated within 
LDA+ U with U =4 eV is consistent with the experimen- 
tal value of 5=187 GPa [16J. For the dependence of spin 
moments of Cr ions on U shown in Fig. 2(d), we see 
clear increasing amplitude of magnetic moments with U 
for both FM and AFM phases. Our calculated value of 
1.47 fiB using LDA for AFM phase is in good agreement 
with previous LMTO calculation value of 1.414 [Ib [29| . 
At a typical value t/=4 eV, the LDA+ U and GGA+ U 
give magnetic moments of 2.46 and 2.62 /j,b for AFM 
PCO, respectively, both of which exceed recent experi- 
mental value of 1.70 [1b 15]. Similar trend has also been 
exhibited in study of BiFeOs [f|. Therefore, inclusion 
of on-site Coulomb energy for adequately describing the 
structural and magnetic properties is crucial. In study of 
Cr02, x-ray absorption and resonant photoemission spec- 
troscopy support the importance of Coulomb correlations 
30]. Huang et al. concluded that the on-site Coulomb 
interaction energy of Cr02 is 3-4 eV through comparing 
their experimental measurements and LDA+ U calcula- 
tions. They found that the shift of Cr 3d spin- up DOS 
slightly away from the Fermi level increases the Cr spin 
moment. In our present study, similar shift of Cr 3d DOS 
is also observed (see below). Overall, comparing with the 
experimental data, the accuracy of our atomic-structure 
prediction for AFM cubic PCO is quite satisfactory by 
tuning the effective Hubbard parameter U in a range of 
3-4 eV within the LDA+ U approach, which supplies the 
safeguard for our following study of electronic structure 
and pressure behaviors of PCO. In following study, we 
present our results within LDA and LDA+ U with con- 
stant [7=4 eV. 



B. Phase transition analysis 

After previous analysis, one question arises: What's 
the ground state of experimentally observed "PbCrOa" 
at ambient condition. In this section, we will test care- 
fully the ground state of "PbCrOs" by calculating the 
total energies of PCO, PCO-CPO, and CPO in cubic 
perovskite structure and some possible distorted per- 
ovskite structures. Their relative energies in G-AFM con- 
figuration calculated within LDA+t/ are shown in Fig. 
3. Clearly, although the equilibrium lattice parameter 
(3.998 A) for PmZm phase of PCO-CPO accords well 
with experimental data of Phase I [l6[ , the Pm3m phase 
of PCO-CPO and CPO are not energetically favorable in 
the whole range of volume. After testing various possible 
crystal structures and atomic arrangements, we find that 
R3c phase of PCO and i?3 phase of PCO-CPO possess 
same energy level with that of cubic PCO in some vol- 
umes. In the whole range of volume, i?3c PCO is more 
energetically favorable than cubic PCO. Considering this 
material only can be synthesized under high temperature 
and high pressure, the high-pressure phase of PCO can 
be understood as crystallizing in Pm3m or R3c struc- 
ture. Increasing the cell volume, as shown in Fig. 3, 
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the R3 PCO-CPO becomes more energetically favorable 
than the two high- pressure phases. This clearly indicates 
that the structure of "PbCrOa" at ambient pressure is R3 
PCO-CPO. Therefore, upon compression there will occur 
structural phase transition among the three phases. As 
shown in the inset of Fig. 3, phase transitions of R3 
PCO-CPO to Pmlm PCO at 1.5 GPa and R3 PCO- 
CPO to R3c PCO at -6.7 GPa can be predicted by the 
slopes of the common tangent rule. The former value co- 
incides well with recent experimentally observed value of 
1.6 GPa, while the latter value is smaller than that. As 
a result, the crystal structures of PCO under high pres- 
sure need more works to clarify. Our present study give 
two most possible structures: Pm3m or R3c. Concerning 
the energetics of the transitions from R3 PCO-CPO to 
Pm3m PCO at 1.5 GPa and from R3 PCO-CPO to R3c 
PCO at -6.7 GPa, Fig. 3 clearly shows that per formula 
unit of PCO-CPO needs 97 meV and 458 meV of energy, 
respectively. 

In Table I, we present the calculated lattice parameters 
for PCO in R3c phase and PCO-CPO in R3 structure. 
It can be found that the equilibrium volume of R3c PCO 
is consistent with that of experimental Phase II and R3 
PCO-CPO comparable to Phase I QjJ. Although Xiao et 
al. reported that both the Phase I and Phase II crystal- 
lize in cubic perovskite structure, our calculations are 
different from their observations. Our present results 
need more experimental works to test. For R3c PCO 
and R3 PCO-CPO, our calculated values of -Efm— -E-afm 
(per formula unit) are —80 meV and —62 meV within 
the LDA+ U formalism, respectively, the corresponding 
bulk moduli 5=164 (#=7.9) and 153 GPa (#=4.3) for 
AFM phase, respectively, 5=157 (#=5.6) and 150 GPa 
(#=4.4) for FM phase, respectively, and spin moments of 
Cr ions are 2.79 and 2.92 hb for AFM configuration, re- 
spectively. Averaged to every Cr 4+ ion, the Efm~ Eafm 
of R3 PCO-CPO is about —31 meV, which is almost con- 
sistent with recent experimental report [15j of the AFM 
ordering at Tn^245 K. For the bulk moduli, values of R,3c 
PCO arc slightly smaller than the experimental value of 
Phase II 5=187 GPa 0, while values of R3 PCO-CPO 
are prominently larger than the experimental value of 
Phase I. The abnormal high compressibility of Phase I is 
due to the fact that the Cr0 6 / 2 to Pb0 6 / 2 transition has 
been compressed to occur under low pressure (see below). 
For spin moments, results of FM phase are almost equal 
to the AFM phase for PCO-CPO in R3 phase and PCO 
in both Pm3m and R3c phases. 



C. Pressure behaviors 

The equation of states of AFM P_CO-CPO in Pm3m 
and R3 phases, AFM PCO in Pm3m and R3c phases, 
and the experimental measured pressure-volume data 
from Ref. [fl6|] are presented in Fig. 4. For Pm3m and 
R3c phases of PCO (R3 phase of PCO-CPO), the relative 
smaller volumes calculated in our scheme compared with 



TABLE I: Calculated lattice constant a, rhombohedral angle 
a, volume V, and Wyckoff parameters for PCO in R3c phase 
and PCO-CPO in R3 structure. For R3c phase, the Wyckoff 
positions 2a (x, x, x) and 66 (x, y, z) refer to the cations and 
anions, respectively. In case of the R3 structure, the corre- 
sponding Wyckoff labels are la (x,x,x) and 36 (x,y,z). 
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experimental Phase II (Phase I) originates from the typ- 
ical overbinding character of LDA. From Fig. 4, one can 
find that our calculated volume collapses of Pm3m PCO- 
CPO to Pm3m PCO and R3 PCO-CPO to i?3c PCO at 
experimental phase transition pressure 1.6 GPa is about 
12.4% and 8.0%, respectively. The former value is larger 
than the measured value (9.8%) in recent experiments 
[ill, while the latter value is smaller than that. Under- 
estimation of the volume collapse value, from R3 PCO- 
CPO to i?3c PCO, can be attributed to the experimen- 
tal fact that the CrOg/2 to PbOg/2 transition has been 
compressed to occur under low pressure of around 0.1- 
1.6 GPa in the experimental compound "Phase I". This 
kind of partial transition leads to abnormal high com- 
pressibility of Phase I compared with CaCr03, SrCr03, 
and high-pressure Phase II of PbCrOs [H, [l?! ■ 



D. Electronic structure 

Figure 5 shows the total density of states (DOS) as well 
as the projected DOS for the Cr 3d, Cr 4s, and O 2p or- 
bitals for AFM PCO in Pm3m and i?3c phases and AFM 
PCO-CPO in R3 phase at selective values of U within 
LDA+ U formalism. Corresponding band-structures cal- 
culated with U—4: eV are presented in Fig. 6, where 
both spin-up and spin-down results are plotted. Since 
spin-down results for PCO in Pm3m and R3c phases 
are same with their spin-up results, as indicated in Figs. 
6(a)-6(b), we only plot in Figs. 5(a)-5(d) the spin-up 
results. For AFM PCO-CPO in R3 phase, slight dif- 
ferences between spin-up and spin-down can be seen in 
Figs. 5(e) and 6(c). Overall, results of PCO in both 
Pm3m and R3c phases indicate that the AFM PCO is 
metallic without accounting for or after switching on the 
on-site Coulomb repulsion [see Figs. 5(a)-5(d) and 6(a)- 
6(b)]. This fact conflicts with the experimental obser- 
vations that the "PbCrOa" is AFM semiconductor with 
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0.27 eV or 0.11 eV activation energy in different temper- 
ature ranges [j| [l5| . In our calculations even increasing 
the amplitude of U up to 8 eV, the metallic state has 
not changed for these two phases of PCO. The metal- 
lic ground states have also been observed for NM and 
FM phases. Besides, inclusion of the spin-orbit coupling 
(SOC) and noncollinearity also can not open a gap at 
the Fermi level. Thus, we only can conclude that the 
high-pressure phase of PbCrC>3 is a conductor. Addition- 
ally, results of CPO and PCO-CPO in Pm3m phase also 
show that they are conductors. Although a gap is opened 
with the Hubbard U=6 eV for Pmlm PCO-CPO (not 
shown), the calculated insulating band gap (2.12 eV) is 
prominently larger than the experimental values 0, [H[ . 
In present study, we prefer to believe that the LDA+ U 
with U =4 eV can give a correct depictions of the ground 
state electronic structures for PCO and PCO-CPO. For 
AFM PCO-CPO in R3 phase, our LDA+ U with [7=4 eV 
calculation open an insulating band gap of about 0.48 eV 
[see Fig. 5(e) and 6(c)]. Figure 6(c) clearly indicates that 
the valence band maximum (VBM) appear at Z point and 
conduction band minimum (CBM) at T point in BZ. We 
find that both VBM and CBM have predominant O 2p 
state character mixed with significant Cr 3d contribu- 
tion. Although the calculated value of band gap (0.48 
eV) is almost two times of the experimental value 0.27 
eV 0, UH , since the activation energy of the Phase I per- 
ovskite increases with lowing temperature , we believe 
that our calculation (valid only at K) can give a proper 
depictions of the electronic structures for low-pressure 
phase of PbCr03. Our calculations clearly illustrate that 
the "PbCrOa" will occur an insulator-metal transition 
together with the phase transition of R3 PCO-CPO to 
Pm3m/R3c PCO upon compression. 

We note that the conductivity of SrCrOa and CaCrOa 
exists controversial [l7l - [2l| . Theoretical calculation for 
CaCrOs [2(| predicted metallic ground state with LDA 
and insulating state with LDA+U. For SrCr03 [2l|, the 
metallic ground sate was observed either with LDA or 
LDA+ U. Besides, using LDA+ U method with U =4 eV, 
Lee et al. [2l| successfully predicted an orbital-ordering 
transition from i§ 9 to d^yldxzdyz) 1 for SrCrOs. No evi- 
dence of orbital ordering within the t^g shell for CaCrOs 
was observed [l9|. In our study of PCO, CPO, and PCO- 
CPO compounds, this kind of orbital ordering in Cr i 2 g 
states has also not been found. 

As shown in Fig. 5, inclusion of the on-site Coulomb 
repulsion will lower the occupation energy of Cr 3d elec- 
trons from the Fermi level and elevate Cr 3d, Cr 4s, and O 
2p orbitals occupation levels near —7.0 eV to high level. 
As a result, localization pictures of electrons occupation 
appear after introducing the Coulomb repulsion. From 
Figs. 5(c)-5(e), one can find that the Cr 4s contribution 
is limited. When Cr ions combining with O ions to form 
covalent/ionic bonds, part of Cr 3d and Cr 4s electrons 
will transfer to O 2p orbital. This kind of electron trans- 
fer behavior can be read from the partial DOS pictures. 
Considering the effect of d-hole creation due to 3d-4s hy- 
bridization, we have also examined the effect of Coulomb 



repulsion on Cr 4s-shell. However, no difference of Cr 
4s orbital has been found for introducing the Coulomb 
repulsion into the d-shell or the s-shell. In Figs. 5 and 6, 
only d-shell is considered to participate in the Coulomb 
exchange. In the whole energy domain, electronic struc- 
tures of AFM PCO in Pm3m and R3c phases have no 
evident differences. The main occupation at the Fermi 
level is from Cr 3d and O 2p orbitals. Our results calcu- 
lated at U =0 eV are consistent with previous calculations 
[29l | . A clear hybridization of Cr 3d and O 2p orbitals 
in the energy range from —7.3 to 0.3 eV can be observed 
at U =0 eV. After switching on the U to 4 eV, this hy- 
bridization energy range is moved to from —6.2 to 0.2 
eV. A well resolved peak of Cr 3d state at around —0.3 
eV at U =0 eV is flatted when the Habburd U param- 
eter being increased to about 4 eV. In addition, a band 
gap in the conduction band is apparent under U =0 to 6 
eV. This band gap increases from 0.6 eV at U =0 eV to 
about 2.0 eV at U—4 eV. The main occupation in the 
conduction band is from Cr 3d orbital with some con- 
tribution from O 2p states. For AFM PCO-CPO in R3 
phase, hybridization of Cr 3d and O 2p orbitals in the 
energy range from —6.3 to —1.0 eV is clear. One narrow 
peak locates just below the Fermi level. 



IV. CONCLUSIONS 

In conclusion, the ground state properties as well as the 
high pressure behaviors of PCO, CPO, and PCO-CPO 
compounds were studied by means of the first-principles 
DFT+ U method. By choosing the Hubbard U param- 
eter around 4 eV within the LDA+ U approach, FM 
and/or AFM ground states were achieved and our calcu- 
lated volumes, bulk moduli, spin moments, and equation 
of states are in good agreement with recent experiments. 
While the PCO-CPO in R3 phase is consistent with the 
experimental low-pressure Phase I, both Pm3m and i?3c 
phases of PCO coincide well with high-pressure Phase 
II. Specially, the semiconductor nature of R3 PCO-CPO 
is in good agreement with experiments. These obser- 
vations explicitly indicate the existence of strongly cor- 
related electronic behaviors in these compounds. Our 
electronic spectrums illustrate a clear hybridization of 
Cr 3d and O 2p orbitals in wide energy range. In con- 
trast to SrCrOs, the orbital-ordering transition from t\ g 
to dl :y (d xz dy Z ) 1 has not been found in these materials. 
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Figure captions: 

Fig. 1: (Color online) Pictorial illustrations of (a) cu- 
bic PCO and (b) cubic PCO-CPO in AFM configuration. 
For PCO-CPO, the first atom along the [111] diagonal 
direction is labeled as Pbl, second Crl, third Cr2, and 
fourth Pb2. One can find that along [111] direction in 
PCO the Pb and Cr alternate while in PCO-CPO the 
Pb and Cr alternate in pairs. Compared with PCO, in 
PCO-CPO half percent of the A site and B site atoms 
are exchanged while in CPO the A site and B site atoms 
are totally exchanged. 

Fig. 2: (Color online) Dependences of (a) total-energy 
differences (per formula unit) between FM and AFM 
(Efm—Eafm), (b) lattice parameter, (c) bulk modulus, 
and (d) spin moments of Cr ions on U for AFM PCO in 
Pm3m structure. 

Fig. 3: (Color online) Comparison of relative energies 
of two unit cells of AFM PCO in Pm3m and i?3c phases, 
one formula unit of AFM PCO-CPO in Pmlm and i?3 
phases, and two unit cells of AFM CPO in Pm3m phase 
vs the volume. All results are calculated within LDA+ U 
at [7=4 eV. Phase transitions of R3 PCO-CPO to Pmlm 



PCO at 1.5 GPa and i?3 PCO-CPO to i?3c PCO at - 
6.7 GPa can be predicted by the slopes of the common 
tangent rule, as shown in the inset. 

Fig. 4: (Color online) The P-V relations of the AFM 
PCO-CPO in Pm3m and i?3 phases as well as AFM PCO 
in Pm3m and R3c phases computed in the LDA+ U for- 
malism. Experimental results from Ref. 161 are also 
presented. The volume collapses at experimental phase 
transition pressure 1.6 GPa are labeled. 

Fig. 5: (Color online) The total DOS for AFM PCO 
in Pmim and i?3c phases as well as AFM PCO-CPO 
in i?3 phase computed in the LDA+ U formalism with 
selective values of U. The projected DOSs for the Cr 3d, 
Cr 4s, and O 2p orbitals are also shown. In panel (e), 
both spin-up and spin-down results are presented. The 
Fermi energy level is set at zero. 

Fig. 6: (Color online) Band-structures of AFM PCO 
in Praira and i?3c phases as well as AFM PCO-CPO 
in i?3 phase computed in the LDA+ U formalism with 
U =4 cV. While the solid lines show the spin-up results, 
the dashed lines stand for spin-down. The Fermi energy 
level is set at zero as shown by the short-dashed lines. 
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